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Advection of a passive scalar ^ in d = 2 by a large-scale velocity field rapidly changing in time 
is considered. The Gaussian feature of the passive scalar statistics in the convective interval was 
discovered in Here we examine deviations from the Gaussianity: we obtain analytically the 
" " ' simultaneous fourth-order correlation function of 6. Explicit expressions for fourth-order objects, 

^ ; like {(6li - Oif) are derived. 

^ . PACS numbers 47.10.-fg, 47.27.-i, 05.40.-fj 

< ■ 

Advection of a passive scalar by an incompressible turbulent flow is one of the classical problems in the theory of 
, turbulence. The problem is related to statistics of temperature or impurities in the flow. The dynamics of the passive 
scalar is governed by 

^! (9t + u„Va - kA)0 = , (1) 

o : 

where the velocity u and the pumping (p are random functions of t and r, k is the diffusion coefflcient. Correlation 
00 functions of 9 should be treated as averages both over statistics of (j) and u. Batchelor [|j was the first who considered 
] the problem for a long-range velocity field. He found the pair correlator of the passive scalar in the case where velocity 
field is very slow. Kraichnan ^ considered the pair correlator in the opposite case of a velocity field changing in time 
very rapidly. A theory for any finite correlation time of the velocity field was proposed in Q . It was proved there that 
Ch ' whatever be the statistics of the velocity field, the statistics of the passive scalar in the convective interval approaches 
Gaussianity as one increases the Peclet number (the ratio of the pumping scale to the diffusion one) . 

In the present letter we aim at finding explicitly the fourth-order correlation function of the passive scalar in 
O ' Kraichnan regime. We assume that the source is (5-correlatcd in time and spatially correlated on scale L: 

: {4>{h,ri)cp{t2,Y2))^5{t^-t2)x2{ri2), (2) 
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where X2{t) tends to zero at r ^ L. We choose the simplest smooth form of the pumping 

X2{r)^P2LV{L^ + r^). (3) 

where P2 is the production rate oi 9^ . A variation of the shape of £,2{r) at r > i will keep all the values of interest 
to be intact at r L. The statistics of a velocity (5-correlated in time is completely defined by the pair correlation 
function which at dimension d — 2 equals 

{Uc{ti,ri)up{t2,r2)) = S{ti - 12) [DLISc/3 - JCa/iiru)] , (4) 
/Ca/3(r) - D{35,,fjry2 - r^rp) , (5) 

where the incompressibility condition V • u = is taken into account. Here _D is a the characteristic strain describing 
the strength of the velocity field, L„ is the velocity correlation length (the size of the largest vortex) which is assumed 
to be the largest scale in the problem. The expression (||) is correct at r -C i„. 

Simultaneous correlation functions of 9 satisfy separate linear equations [§j^J^ . The equation for the pair correlation 
function /(|ri — r2|) = (^1^2) can be solved explicitly for arbitrary r. For separations r ^ r^, where — is 
the so-called diffusive length, one gets 

/(r) = g{ln(l + LVr2) + (iVr2)ln(l + rVi2)| . (6) 
The equation for the fourth-order simultaneous correlation function of the passive scalar !F = (6*1^2^3^4) is 
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-CT{ri,r2,r3,r4,) = *(ri2,r34) + *(ri3,r24) 

4 

'^{r+,r-) = .f{r+)x2{r-) + f{r-)x2{r+) , 



(7) 

(8) 
(9) 



where = |r. 



arbitrary dimension d 



The equation (|7|) can be generahzed for an arbitrary multi-scale velocity field and for an 
HJ^. Although at any d r.h.s. of decomposes into three parts each depending only on two 
separations, only in the case of the large-scale velocity field the operator £ supports the decomposed form: acting 
on an arbitrary function of the two vectors it produces a function of the two vectors too. The corollary gives us a 
possibility to establish the following general form of the solution 



Jp(ri, r2, ra, = i^(ri2, r34) + F(ri3, r24) -I- -F(ri4, r23) 



(10) 



Since we are looking for a solution in the restricted region r <^ Lu one could add to (0) a zero mode of the operator 
C. However, the equation (Q) stems from the dynamical approach (see |^for details) that keeps no room for another 
solution besides of the decomposed one. The function i^(r_|_,r_) from (1^) satisfies the following equation 

-[£' + £d]F(r+,r_) = <i>(r+,r_), 



(11) 
(12) 



In ( p^ we passed to the variables ^ = ln[L^/(r+r_)], y = ln[r_/r+] and -d = arccos[(r+r_)/(r+r_)]. The physical 
boundary conditions imposed on F are: _F'(r_|_,r_) — » at r+ or r_ oo. 

As was shown in |l], correlators of 9 do not feel diffusion if r± ^ r^. Thus at r we can omit the diffusive 

term in ( pi] ) that together with the independence of <&(r_|_, r_) on allows us to write out a solution of ( pl| ) as a sum 
F — F+ + F_ , where F± satisfy the equations 



-2Dsm'^[dl + d^F±=fir±)X2ir^) 



(13) 



Here r± = \/S/\ sinz9| exp(=Fj//2) and 5", z?, y should be treated as independent parameters. We conclude that, by 
construction, F+{S, "d, y) = F-{S, i), —y). Taking into account symmetry properties of (^3|) we find that F± is invariant 
under d —i) and under — > tt — i?. To solve the equation (|l3|) one can use the resolvent TZ of the Laplacian + dy 
figuring in ([l^): 



- + 92] ^(^^ ^',y-y')^ S{{) - d')S{y - y') , 



(14) 



One should impose the zero boundary conditions on the resolvent at the boundaries of the strip: < < tt and 



—00 < y < +CXD, since in accordance with the definition r± 
written explicitly: 



if "i? — > 0, TT or y ■— > ±cx). The resolvent can be 



7^ = — In 

An 



sinh^(y/2 - y'/2) + siTi^{'d/2 + /2) 



sinh^(j//2 - y'/2) + sin^(i9/2 - ir/2)_ 



Convolution of the r.h.s. of (O) with the resolvent gives 



F+ = 



4Dr^_ sin^ d 



dwf{rj^w)T ^ In 



{w+ + tY + w'^ cot^ 1? 



(15) 



(16) 



where w± = \w ±w | , t 
F+ a.td = 0: 



L'^rJ^ sin -d + w ^. Note that this formula leads to the weak angular singularity of 



d^F+\^^^ - -KPi/D^e^y ln(l + e"^) - e^) ^ . 



(17) 



Below we will show that the diffusion smoothes the singularity at the smallest angles. Generally, ( ^q ) substituted into 
F = + Fj^{y —y) and further into ( p^ ) closes the problem of finding the fourth-order correlation function T of 
the passive scalar for all separations being larger r^. 

At L ^ r± the integration in the formula dlq) can be performed, then 
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p2 

F = j^{\n'[L/r+]/2 + \n'[L/r^]/2 + \n[L/r+] + ln[L/r_] + 

+y{z)+y{z) + Oir±/L), (18) 
yiz) = ^{7rlm[e--] -Re[l7(-ie-)] -Re[e-2-f](ie--)] , 

n{u) = — +ln^u + Li2(l - m) + Li2(l + m)} , (19) 

where z — -d + iy and z = i9 — ij/, and Li in ( p^ ) is the integral logarithm. Let us emphasize, that the second part of 
(|l|), oc [y{z) + y{z)], turns out to be a zero mode of the Laplacian since it is the real part of an analytic function 
of z. The general result (^) shows a remarkable coincidence of F within a logarithmic accuracy with its coUinear 
limit F{'d = 0). It follows from the dynamical approach described in ||l[. It was demonstrated there that F is an 
average over velocity's statistics of the product of times required for the vectors r^- and r_ respectively to grow upto 
L if their edge points move along Lagrangian trajectories. An evolution of the vectors, being originally deep inside 
the convective interval r± <C L, can be divided into two different stages: the first one is a quick collinearization of 
the geometry, and second one is a long coUinear stretching giving the major contribution into the correlator. At 
L » r+ > r_ (|l|) gives 

[F - /(r+)/(r_)] « ^{ln[L/r+] + I + - 7r|^|/2} . (20) 

Here in the l.h.s. we subtracted from F the major Gaussian part. 

Let us now return to F^ to examine its behavior at small angles. The angular singularity (^^ is formed by the 
z-dependent part of ([l6| ) that is a zero mode of £' . At small angles it should be replaced by a zero mode Z of C' + Cd, 
where the diffusive operator Cd can be rewritten in terms of the set of variables (■&, y, ^) as 

Cd = 2k cosh[y] {9| + 8^ + + 4 tanh[2/]a^9j, } . (21) 

The angular singularity is formed by the z-dependent part of (|l6|) which can be treated as a zero mode of £' by 
analogy with (|l^). At small angles it should be replaced by a zero mode Z of the sum C' + Cd- We are looking for a 
zero mode Z that is proportional to |??| at values of the angle, where it is possible to neglect the effect of diffusion at 
all, and is regular at i9 — > 0. At the smallest value of i9 d^Z is much larger than dyZ and d^Z. Thus, Z should satisfy 

{ (13252 + 2^d^d^ -d^+dl)+Ae^dl}Z = 0, (22) 

where A — KCOsh{y)/ (16DL^) can be considered as an arbitrary parameter. Taking into account ( p7[ ) we find the 
following solution of ( p^ possessing the desirable behavior 

Z « TT-^ [e^y ln(l + e-y) - e^] sjd^ + k coshye^ {DL^) . (23) 

The characteristic angle -iJo = \/ kD~^L~^ cos]i[y] exp(^/2) turns out to be small at r± ^ rd- Using the explicit form 
(^) we conclude that for 1} ^ do: d^Z ^ d^^Z, dyZ ~ Z, d^Z ^ Z. Those estimations justify above calculations. 
To summarize, at r± 3> rd the diffusion is relevant only at small angles ^ i?o where it influences angular derivatives 
of F but it gives a negligible correction to the expression (|l^). 

As long as one considers a correlation function at sufficiently small distances, say P at r„ < the account of 
diffusion is unavoidable. The major value of the function can be found directly from the suitable expression in the 
convective interval by putting there, in terms formally divergent at r — > 0, rd instead of the smallest distances (see 
Appendix A of for the proof). Thus, using (|l^), we find with the logarithmic accuracy the following fourth-order 
objects 

at L > ri3,ri4,r23,r24 > > > Td, ri2||r34 

{(^1 - O^Y) « 12^(ln2[ri2/rrf] +ln[ri2/rd]) , (24) 

m - e^m) {{9, e,f){el) « 2^ Hr.^/rd] , (25) 

((^1 - ^2)'(^3 - ^4)') - ((^1 - e2f)m - Oa?) ~ 4^ ln[r34/rd] . (26) 
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(gj shows that with the logarithmic accuracy (/(O) « f{rd)) ((6'i-6'2)^) = 2/(0)-2/(ri2) « 2P2/ Dlii[ri2/rd]. Thus, 
it follows from (|2^ ) that the flatness ((6*1 — 6'2)'^ ) / (((^ 1 — (^2)"^)) ~ 3 in accordance with the asymptotic Gaussianity 
proved in To overcome from the correlators ( [25|j2^ ) to respective correlators of the dissipative field e — k{V6)'^ we 
should differentiate them over ri2,r34 and replace the separations after all by r^. One gets the following estimations 

at L»ri3»rd {{ei9l)) = K{{{Vei)^9l)) ^ / D , (27) 
((6i63))=^2(((V0if(V03n)-P|- (28) 



Unknown multipliers behind the parametric dependencies in the formulas ( |27|j2§| ) are of the order of unity. We 
postpone an explicit calculation of the multipliers, that requires a direct account of diffusion. The formulas (27|2^) 
show (in accordance with M) the zero dimensionality both for the passive scalar ^ and the dissipation field 

Note that in the region Lu ^ rij,ru ^ L,rik,rii,rjk,'rji the found correlator T does not decay. The decay occurs 
only at the largest scales ri2 , ^34 ^ L„ which are out of the scope of the present study. 

The proposed solution for the fourth-order correlation function of 9 could serve a starting point for studying passive 
scalar correlations in the general case of a multi-scale (but short-correlated in time) velocity field, the problem is 
introduced by Kraichnan g. A closure applied recently to the model, yields an anomalous scaling, particularly, 
for ((^1 —6*2)"') that conflicts with an analytic consideration 0|. Future attempts of a perturbative study based on the 
results of the present letter could solve the collision. From another hand, an inclusion of the present results into the 
scheme proposed in |10 1^ could help in better understanding the problem of the direct cascade in two-dimensional 



turbulence. 
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